A simple vision of current induced spin torque in domain walls 
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The effective spin pressure induced by an electric current on a domain wall in a ferromagnet is 
determined using a simple classical model, which allows us to extend previous theories to arbitrary 
domain wall widths. In particular, the role of spatially non-uniform components of the torques are 
analysed in detail. We find that the effect of the current is mainly to distort the domain wall which 
should enhance de-pinning. We also find that in the limit of thin domain walls the current-induced 
pressure changes sign. 



Traditionally, spin electronics has been dealing with 
resistance changes induced by different magnetic con- 
figurations, like those encountered in 'Giant Magneto- 
Resistance' and 'Tunneling Magneto-Resistance' devices. 
This has led to important developments in information 
technology applied in computer read-heads and the more 
recent magnetic memories (or MRAMs) for high den- 
sity information storage. For the latter, one has to be 
able to read and write the magnetic information. Writ- 
ing is presently achieved by local application of a mag- 
netic field using current lines crossing on each individual 
cell. This procedure is the main source for cross-talk be- 
tween magnetic elements because of field leakages which 
can influence neighbouring memories. Thus, it would 
be most useful to switch the magnetic configurations di- 
rectly with a current flowing in the spin-valve. The rel- 
evant effect is called spin-torque and it has been pre- 
dicted by Slonczewskii and Berges^ and observed in the 
late 90—. It is now well-known that magnetization can 
be reversibly switched in ferromagnetic-normal metal- 
ferromagnetic trilayers by a large current crossing the in- 
terfaces. Besides, it appears that this effect can also pro- 
duce microwave-frequency oscillations of the thin layer—, 
possibly allowing coherent microwave generation by an 
assembly of nanopillars^. An alternative way to control 
magnetic configurations is to move domain walls (DWs) 
with a current as first proposed by Berger— and demon- 
strated a long time agoi Recently, theoretical studies^ 
and numerical calculations^ have shown that a large part 
of the current-induced torque does not push the DWs and 
the pressure only originates from spin-flip events relaxing 
the conduction electrons' spins. 

We propose here a simple picture of the relevant 
physics involved in the generation of the different compo- 
nents of the torque in domain walls. The basis for it lies 
in understanding the evolution of the spin of a conduc- 
tion electron as it crosses a domain wall and its reaction 
on the local moment. The relevant interaction is that 
between localized (the local magnetisation) and delocal- 
ized (conduction) electrons which can be expressed with 
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The exchange interaction splits the conduction electrons 
in two populations with spins parallel (up) or antiparallel 
(down) to the local moments. A current generates a plane 
wave of electrons whose wave functions can be expressed 
as a spinnor with two components, up and down, trav- 
elling with different wave vectors. The question is then 
to analyse how these two components evolve when the 
electrons are forced to cross a region where localized mo- 
ments change direction in space, i.e. a domain wall. The 
proper way to do this is to write Schrodingcr's equation 
and match wave functions and their derivatives at the 
borders of the DW. Solutions of the problem are complex 
for any particular value of the domain wall width, but an- 
alytical solutions can be easily found for the two limits 
of an abrupt, or a very long domain wall (see for exam- 
ple ref^S). There is another way of treating the problem 
which is to consider conduction electrons as free particles 
entering a region in space where a local field changes di- 
rection (the DW). Their spin evolution is then obtained 
by writing the Landau-Lifshitz equation. This is sim- 
pler, but one has to overcome conceptual problems linked 
with the nature of the electrons crossing, which are band 
particles. These two visions of the problem are actu- 
ally equivalent in the limit where the amount of reflected 
wave can be neglected^. Hence, this is valid for a DW 
width much greater than the Fermi wavelengthii. Sev- 
eral works in that field (addressing more particularly DW 
resistance) have chosen either the 'wave' approachi2iiii4 
or the 'particle' on o 15 i 16 . Because in the wave approach 
one has to use perturbation theory, spin torques have 
been calculated in the limits of wide, or very narrow do- 
main walls. Here, we propose to use the simplicity of 
the 'particle' approximation to describe the spin evolu- 
tion in the domain walls and to extend previous results 
to intermediate values of the DW width. 

Classical model: 

If both magnetization and electron magnetic moment 
(p = — q^bs) are considered classical vectors, the spin 
dynamics obeys the basic precession equation: 



the exchange interaction, < S > / S — —MJM S refers to 
localized spins, and s to the conduction electrons spins. 
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where m = M /M s is the unitary vector of magnetization. 
This is similar to the classical Landau-Lifshitz equation 
with the external field replaced by the exchange field, 
which reflects the difference between the s-d Hamiltonian 
considered here and the simple Zeeman one. 




FIG. 1: Bloch domain wall : conventions 

We consider here the frame moving with the elec- 
tron crossing the DW, in which the magnetization varies 
continuously in time. For simplicity, the magnetisation 
rotates in the wall plane (Bloch type) but the follow- 
ing derivation would equally apply to a rotation around 
any axis (as long as the system can be considered one- 
dimensional). The moving frame is defined by the three 
vectors (e^,e"e,ey) as shown in figure^ The evolution of 
the electron magnetic moment is thus simply described 
by: 




FIG. 2: Evolution of the conduction electron magnetic mo- 
ment in the frame rotating with the magnetization in the 
domain wall. The precession of ft, is not around the local mag- 
netization but around the effective field direction m/T ex —6e y 



of the magnetisation and its rotation vector: m/T ex — 9e y . 
Every precession, the electron's magnetic moment lags 
from the magnetization and then catches upi&. As a 
consequence, the average over time of the mistracking 
between fl and the magnetization is not zero as it lies on 
the effective field: 
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Defining r ex = h/SJ ex we get: 
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As the conduction electrons cross the DW, their spins 
precess with the so called Larmor period. It is conve- 
nient to consider first a long DW (where 9r ex 3> 1) with 
a linear variation of the magnetization angle 9 (the sec- 
ond derivative of the angle is neglected 9 = 0, and the 
electrons spin remains mainly aligned with the magneti- 
zation), which leads to the simplified equations: 



This is a crucial point as it explains why, even if we 
consider a large number of Larmor oscillations of fl in the 
wall, the resulting effect on the magnetization will not 
be averaged to zero. Because at any time during DW 
crossing the total magnetic moment has to be preserved, 
a torque is induced on the magnetization by each electron 
(per unit volume): 
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This part is adiabatic as there is no energy loss in the 
process. The resulting torque can be decomposed into a 
constant and a periodic terrain. For long DWs, the pe- 
riodic part averages to zero and the constant spin torque 
term: 




1 

He + -5- Ate 

T 

' ex 
1 

Ms/ + ZT^y 



9 guB 



(4) 



In fact the magnetic moment fl is precessing around an 
effective field (fig-EJ whose direction is given by the sum 
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is not in a proper direction to move the DW^. Re- 
cently, it was shown that the term responsible for current 
induced DW motion is due to a non-adiabatic part of the 
torquo^. In order to account for it in our classical model, 
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it is necessary to include the spin-flip relaxation time: 
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where fl eq = g/iB/2e r is aligned with the magnetiza- 
tion. In usual ferromagnetic materials, the spin-flip re- 
laxation time T s f is large compared to the Larmor period 
Tl = 2ttt £X , and the effect can be considered a pertur- 
bation. Then the average mistracking of the electron 
magnetic moment fj, y as derived in cq. [S] can be kept. 
The induced effect is given by the reaction torque of this 
spin-flip related term on the local magnetization: 
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Under an electric current density j = je y , and choosing 
to introduce the magnetization gradient in the wall we 
then get for the two torques: 
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Where the polarization P is added as a prefactor to ac- 
count for the partial polarization of the charge carriers. 
So the adiabatic spin torque accounted by various au- 
thorsi2*i£ii2ii*i can be explained as the reaction from the 
precession of the conduction electrons magnetic moments 
around the effective field. The second term originating 
from the spin flip scattering has an amplitude reduced 
by a factor r ex /r s / typically around 1/30. When writing 
the micromagnetic equations of motion of the DW, one 
can see that only this (small) second term applies a pres- 
sure that pushes the wallS*i^. This can also be seen in 
the following manner: from the Landau-Lifshitz equation 
it can be seen that the precession of a conduction elec- 
trons' spin around the magnetization is equivalent to a 
field oriented along the rotation vector acting on the mag- 
netization. Because the main spin evolution of electrons 
crossing the wall is a rotation following the DW mag- 
netization, the equivalent field they generate is directed 
perpendicular to the plane of the DW and hence to the 
magnetization of the domains. As a result, this does not 
push the wall in any specific direction but only induces 
a canting of the wall magnetization. On the other hand, 
spin-flip processes tend to align the conduction electrons 
spins with the local magnetization. The effective field 
this generates is perpendicular to both the local magne- 
tization and the initial direction of the electron spin as 
the event occurs. Considering average quantities during 
DW crossing, the local magnetization lies globally along 
its direction in the middle of the wall and the electron 
spin lags behind along the direction of rotation of M (that 



given by 9). Thus, on average, this effective field is per- 
pendicular to both vectors, i.e. it is along the direction of 
the domains magnetization. Therefore, this is applying 
a pressure which tends to push the DW. 

Hence, our classical model leads to deformation and 
pressure expressions which are consistent with the results 
of semi-classical theories. The advantage of our simple 
formalism is twofold. It is now possible to study in detail 
the spatial evolution of the torque along the width of the 
DW, and it also makes possible to explore the torque 
for wall widths of the order of the precession length: 
A ~ \l. Indeed, the only essential approximation of 
the model is that the reflected part of the plane elec- 
tron wave impinging on the DW can be neglected. This 
imposes DW widths much larger than the Fermi wave- 
length of the electronsiSiii, which is not a very stringent 
condition since the latter is around 3 A in metals. 

Numerical calculations for thin walls: 

We have numerically calculated the current-induced 
torques exerted locally on the wall from eq. [SI without the 
simplifications associated to large domain walls. Pressure 
and deformation terms are labelled respectively T p and 
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These are also averaged on the different directions the 
Fermi velocity can take on the Fermi sphere. This in- 
duces decoherence of the torques as electrons travelling 
with different components of their velocity on the direc- 
tion perpendicular to the DW, v y , have different Lar- 
mor precession lengths. The resulting effect is a damp- 
ing of the averaged spin precession (hence the oscillat- 
ing part of both torques) after a few Larmor periods 
2irT ex . This allows us to extend the work of Waintal and 
ViretiS who derived the periodic torque but did not in- 
clude the spin-flip scattering, which lead to a zero average 
of T p . Figure |21 presents numerical computations of the 
torques along different walls: linear and Bloch type, with 
A = 5A L or A = X L (with X L = and A s/ = 50\ ex . 

An obvious result is that the exact DW shape is impor- 
tant for the torques as also pointed out by Xiao et aliSS. 
For a linear wall (dashed line) where 9 is discontinuous, 
oscillations are enhanced compared to those in a Bloch 
wall (solid line). Thus, in the pure Bloch walls of bulk 
materials, the periodic torque is small but in short lin- 
ear walls, like these expected in constrictions^, periodic 
components are significant. In real systems, where DWs 
are often pinned on defects or impurities, large torque 
oscillations can also be expected because of the abrupt 
perturbations defects have on the local magnetization. 
This is schematically shown in fig. 0] for a Bloch wall 
pinned on a non-magnetic impurity. Beyond the simplic- 
ity of the chosen system, the result of these large torque 
oscillations should be an efficient depinning of the walls. 
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FIG. 3: "Distortion" and "Pressure" torque parameters Td 
and r p in the domain wall under a current density j with 
A = 5Al (top curve) and A — \l (bottom) and = 1/30. 



Bloch walls. This reflects the different relevant length- 
scales associated to each wall: in a linear wall, because of 
the abrupt border of the DW, the Larmor precession is 
maximum and it becomes the only relevant lengthscale 
whereas for a smooth Bloch wall, the periodic part of 
the torque is smeared and the wall width becomes more 
relevant. Hence, when the DW width is larger than the 
spin diffusion length, the torque decreases because the re- 
laxation reduces the average mistracking, and hence the 
pressure torque. As a result, the asyptotic values for v p 
in fig.[5]are — 4 s2 - for linear walls and zero for Bloch walls. 
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FIG. 4: Numerical simulation of the effect of a non-magnetic 
impurity on the torque v p , for a Bloch DW of width A = 
2.5Al. The figure is a superposition of ID simulations for 
the varying width of the non-magnetic material, considering 
it does not affect the domain wall profile. 



Let us now turn to the average torque values on the 
domain wall width A. To allow for a straightforward 
comparison with ref^ we introduce the two velocity-like 
parameters Vd and v p corresponding to the torques ex- 
erted on the magnetization, respectively in the eg and e y 
directions: 
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The dependence of the relevant quantities Vd and v p 
with A is presented in fig. [S] Asymptotic values for 
large domain walls (A 3> Al) are different for linear and 




Domain wall width A (XJ 

FIG. 5: Averaged distortion (vd) and pressure (v p ) as a func- 
tion of the domain wall width A for linear and Bloch walls. 
The component applying pressure(iip) changes sign and in- 
creases steeply for thin walls (the asymptotic value for wide 
linear walls is represented by the dashed line). 

For thin domain walls with width close to the Larmor 
length Xl, the oscillatory torque neither tends to zero 
nor is it washed out by the average on the Fermi sphere. 
Interestingly, as A is reduced, the v p term changes sign 
and becomes much larger than its value for wide walls: 
jP£tB_z^_ This is at first sight surprising, but can be 
understood as follows: when the electrons first enter the 
DW, the magnetization rotation generates a tilt in the ef- 
fective field around which the conduction electrons spins 
precess. Hence, their first period takes them away from 
the local magnetization and the corresponding torque is 
opposite to that generated by spin-flip scattering. It is 
also much larger because at first all the electrons are 
precessing in that same way. In very thin DWs, this ef- 
fect dominates. Thus, in appropriate systems, the ef- 
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fective torque could be an order of magnitude larger 
than that in conventional domain walls. This is particu- 
larly relevant for ferromagnetic materials with very large 
anisotropy, where the predicted domain wall width is only 
a few nanometers, which is close to the Larmor preces- 
sion length Xl- One could also imagine making walls 
thin in nanometer-sized constrictions, or increasing Xl 
as in magnetic semiconductors. This might actually ex- 
plain results in GaMnAs, where DWs have been shown to 
move under much lower current densities^ correspond- 
ing to an efficiency of 30% for the spin torque (an order 
of magnitude higher than in NiFe). This would be the 
case in our calculation for A ~ 2Xl, which is reasonable 
in GaMnAs because of a relatively weak exchange and 
large anisotropy. Interestingly, for thin walls, conduction 
electrons are injected in the magnetic domain beyond 
the DW with their spin significantly misaligned. Spin 
flip events will eventually realign the spins thus leaving 
some angular momentum within a spin diffusion length 
from the DW. It is likely that this will result in large 
magnon emission in the domains. This effect is beyond 
the scope of this paper, but one can conjecture that for 



a large current, the magnetization beyond the DW could 
be destabilized. 

In summary, the model we develop here allows us to 
study the evolution of the conduction electrons spins as 
they cross a domain wall. For long walls, the average 
torque is mainly due to a constant term coming from the 
global spin rotation while following the magnetization 
direction. This tends to distort the DW without push- 
ing it. The only contribution applying a pressure is that 
due to spin-flip scattering events which amounts to a few 
percent of the total torque. For thinner walls, the contri- 
bution from the periodic torque, which depends on the 
exact shape of the DW, becomes important. The total 
torque can be significantly enhanced and it is even found 
to change sign for very thin walls. Thus, we predict that 
for domain walls of width close to the Larmor precession 
length, the current induced spin pressure is an order of 
magnitude higher and in a direction opposite to that for 
large DWs. 
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